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■ Abstract. We study the possible singularities of isotropic cosmological models that have 
a varying speed of light as well as a varying gravitational constant. The field equations 

^ , typically reduce to two dimensional systems which are then analyzed both by dynamical systems 

■ techniques in phase space and by applying the method of asymptotic splittings. In the general 
' case we find initially expanding closed models which recoUapse to a future singularity and open 
, universes that are eternally expanding towards the future. The precise nature of the singularities 

^\ • is also discussed. 
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O^' 1. Introduction 

5-H . Varying speed of light (VSL) cosmologies have recently received considerable attention as 

alternatives to cosmological inflation to provide a different basis for resolving the problems of the 
standard model (see for example and refs. therein). Instead of adopting the inflationary idea 
^ i that the very early universe experienced a period of superluminar expansion, in these universes 

;h ' one assumes that light travelled faster in the early universe. In such cosmological models the 



known puzzles of the standard early universe cosmology are absent at the cost of breaking the 
general covariance of the underlined gravity theory. We do not enter here into a discussion on 
the foundations of VSL theories and the conceptual problems arising from the very meaning of 
varying the speed of light (see P for a discussion). 

In this paper we report on some preliminary results of our on-going work [2j on the application 
of the method of asymptotic splittings introduced in li to study the singularities that may arise 
in VSL cosmologies. In particular, we focus here on the VSL model proposed in |^ and further 
investigated in 0. An important characteristic of this model is that one assumes minimal 
coupling at the level of the Einstein equations so that a time-variable c should not introduce 
changes in the curvature terms in Einstein's equations in the cosmological frame, Einstein's 
equations must still hold. As Barrow points out c changes in the local Lorentzian frames 
associated with the cosmological expansion, a special-relativistic effect, so that the resulting 
theory is not covariant and one has to make a specific choice of time coordinate. Choosing that 
specific time to be comoving proper time the Friedman equations retain their form with c{t) and 
G{t) varying. 



We assume an equation of state of the form p = (7 — 1)/0C^, with < 7 < 2 and we write the 
Friedman equations with varying c (t) and G (t) as follows: 
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Here a is the scale factor, A; = 0, +1, or —1. Differentiating the first and using the second equation 
above we obtain the conservation equation 

p + 3jp- = -p— + 3———, (3) 
a G 47rG 

and setting a = x, x = y we obtain the non- autonomous system 

X = y, 

y = -li(37-2)xp, (4) 

o y A: c{t)c{t) 

P = -37-/5 - P777TT + 3- 



subject to the constraint 



y^ kc^ (t) _ SttG (t) 
x"^ x"^ 3 ' 



subject to the constraint 



2. The fiat case 

For flat [k = 0) models we set z = Gp and the system @ becomes 

X = y, 

Svr , 

y = — — (37 - 2) xz, 



z + 3-f-z = 0, (5) 

X 

y2 ^ Svr 
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This system is exactly that which was analyzed using the method of asymptotic splittings in , 
and we rewrite their result. Putting 

X = ar^, y = Pt'', z = 6t^, 

we find the balance 

2 2 1 

p = —, q = p-l, r = -2, (3= —a, 6 = - — 2 > « = arbitrary. (6) 
37 37 D7r7^ 

Computation of the eigenvalues of the Kovalevskaya matrix K = Df (a) — diag (p) yields the 
values —1,0, Applying the constraint to the second equation we obtain an integrable 

system with general solution 

:, = [At + Cf/^\ y = M(^t + C)2/37-i, , = JL^^^At + C)-\ (7) 

37 D7r7^ 



where A and C are constants of integration. We see that the density function always has a 
finite time singularity. On the other hand, the Hubble parameter blows up in finite time only 
for 7 < 2/3. It is interesting that the leading-order terms © obtained from the calculation of 
the dominant balance already contain all the information provided by the exact solution 0, an 
effect due to the fact that the system is weight-homogeneous. We will report on the details of 
the full analysis of this system elsewhere. 

3. Reduction to two dimensions 

We assume a power-law dependence of the function c, 

c = cqo^, n S M, G = constant, 

and use a system of units with SirG = 1 = Cq. We can avoid having to deal with denominators 
if we set x = 1/a, H = a/a. Then the system (j2|, © becomes 

X = —xH, 

H = -H^-^^^p, (8) 

D 



p = -3jpH + 6knu^-'^''H, 



and the constraint reads 



We use the constraint to eliminate x, and write the system in the form 

6 

p = - {3-/ - 2n) pH - 6nH^ . (10) 

Prom Eq. ® it follows that the phase space of ^ is the set { (if, p) £ : p > 0, p - Si?^ > o} 
for closed models, and {{H,p) ^M? : p > 0, p — 3H^ < O} for the open ones. Note that the 
corresponding equations in general relativity have the feature that the conservation equation 
corresponding to the second equation in is just p = —3jpH, which implies that the line 
p = is invariant, i.e., the trajectories of the system cannot cross the line p = 0. On the other 
hand, here equations (|T?1|) without further assumptions do not guarantee that a solution curve 
starting at a point with p > will not eventually enter the region with p < 0, which of course 
is unphysical. 

A first task of the method of asymptotic splittings is to find all possible asymptotic forms 
H = ar'P, p = 13t^ admitted by the system (|T?H) . One balance gives p = —l,q = — (87 — 2n) 
with a = 1, /3 = arbitrary. This has K-exponents (—1,0). A second interesting balance is for 
p = — l,g = —2 with coefficients given by 

2 ^ 4 \ / 1 6n 



37 V l-n'^ (3^_2)(l-n)' 

The K-exponents are in this case given by the forms 

2 , A / 37 + 2n-2 

-1, — (37 + 2n-2) , -1 



37 / V 



These results lead to formal series expansions of the solutions in a suitable neighborhood of 
the finite time singularity, they will be presented elsewhere, j2j. It is interesting that a phase 



H 



Figure 1. The phase portrait of (|inj) 



space analysis shows that only closed models run into a finite time future singularity. In fact 
the phase portrait of with n = —1/2,7 = 1 is shown in Figure 1. The dashed line in this 
Figure separates the closed from the open models. We see that every initially expanding universe 
starting below this line eventually approaches the origin which corresponds to zero density and 
expansion rate. On the other hand, initially expanding universes starting above the dashed line 
eventually contract and H — oo, /3 ^ cxd in a finite time. Numerical experiments confirm the 
above analysis. 

4. The general case 

The reduction of the three dimensional dynamical system © to the two-dimensional system 
(|1U() is not unique. We can use the constraint © to eliminate p instead of u. However, we shall 
see that assuming a power-law dependence for G as well, we can handle both cases at once. In 
this Section we assume that 



and use a system of units with SvrGo = 1 = Cq. Setting again x = 1/a, H = d/a, the system Q, 
©, © can be written as 



c = 



Cod 



G = Goa 



n, m £ M, 



X 



xH, 



(11) 



H 
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H. 



subject to the constraint 



^2 + A;x2-2« = ^px-"^. (12) 
We use the constraint to ehminate p, and arrive at the two-dimensional system 

X = —xH, 

H = _g2ff2_ (37-2) 2-2n_ 

2 2 ^ ^ 

Note that from (|12l) the phase space of H13|) is the set 

{(x,i/) G : X > 0, + > 0.} (14) 

This system is now in a form suitable for both a dynamical systems analysis and an application 
of the method of asymptotic splittings. Full details will be given elsewhere but we report here 
on some partial results to give a flavor of the analysis. Putting x = ar^, H = fSr^ we find a 
possible balance with q = —l,p = —2/37 and (3 = 2/87 = 1, a = arbitrary, acceptable for n < 0. 
The Kovalevskaya eigenvalues for this balance are at —1,0, compatible with the arbitrariness of 
a. 

In the case where all terms are dominant, the vector field is 



f(0) 



-xH 



and we find the balance (disregarding the case a = 0, /3 = 0,) 



n-lj' 1-n' " (37-2)fc(l -n)2 ' ' 

with K-exponents given by 

2 - 2n - 37 
1 — n 

Note that if n = —1/2 and 7 = 1 (dust), or if n = — 1 and 7 = 4/3 (radiation) the second 
eigenvalue is equal to 0. 

We may proceed to analyze the system ()13p in phase space. First, we write the system it as 
a single differential equation, 

dH _ + (37 - 2) fcx2-2« 

~ 2^ ' ^^^^ 

and we set = z to obtain a linear differential equation for z which is easily integrable. We 
find that has the general solution 

= (37 -2) fc ^2-2n^^^37 if37 + 2n^2, (16) 
2 — 37 — 2n 

^ (37-2)A;x3'^lnx + Cx^^ ff 37 + 2n = 2, 

where C is a constant of integration. Suppose that 37-|-2n 7^ 2. We observe that the leading-order 
terms of the first decomposition, namely 



H 



H 




^ = +1 

Figure 2. The phase portrait of H13() 



reproduce the Cx^"^ part of the solution The leading-order behavior of the last 

decomposition, 

recovers the part of the solution. 

The integral curves of (fT^ allow us to sketch the phase portrait of ((T^ . Let us assume for 
concreteness that n = —l, 7 = 4/3 (radiation). There are two cases to consider: 

• Closed models. The phase space (see CU) is the half-plane x > 0. Equation (fT6|) yields 

= —2x'^+Cx'^, which implies that C > 0. Any orbit starting in the first quadrant satisfies 
X > 2/C, i.e., there are no solutions approaching the origin. For C > 0, any orbit of (|13j) 
starting in the first quadrant crosses the x— axis at 2/C and remains in the fourth quadrant. 
We conclude that any initially expanding closed universe reaches maximum expansion at 
some finite time and then recollapse begins, i.e. H approaches — cxo in a finite time. The 
phase portrait is shown in Figure 2. 

• Open models. From (|14|) we see that the phase space is the set x > 0, > x^. In the 
first quadrant this set is represented as the area above the dashed line in Figure 2 B. Eq. 

yields = +2x'^ + Cx^, which implies that for C > any orbit starting in the first 



quadrant asymptotically approaches the origin. For C < 0, there are homoclinic curves 
connecting the origin with itself. However, we must remember that the allowed initial 
conditions for expanding universes lie above the dashed line in Figure 2. It can be shown 
that these trajectories also asymptotically approach (0, 0) . We conclude that any initially 
expanding open universe remains ever-expanding. 

In all, we find that VSL cosmological models may share many interesting dynamical 
characteristics not fully present in the more conventional models so that a more detailed study 
of the singularities in these universes is worthwhile as a future challenge. 
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